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We study the damage process of fiber bundles in a wedge-shape geometry which ensures a constant 
strain gradient. To obtain the wedge geometry we consider the three-point bending of a bar, which is 
modelled as two rigid blocks glued together by a thin elastic interface. The interface is discretized by 
parallel fibers with random failure thresholds, which get elongated when the bar is bent. Analyzing 
the progressive damage of the system we show that the strain gradient results in a rich spectrum 
of novel behavior of fiber bundles. We find that for weak disorder an interface crack is formed 
as a continuous region of failed fibers. Ahead the crack a process zone develops which proved to 
shrink with increasing deformation making the crack tip sharper as the crack advances. For strong 
disorder, failure of the system occurs as a spatially random sequence of breakings. Damage of the 
fiber bundle proceeds in bursts whose size distribution shows a power law behavior with a crossover 
from an exponent 2.5 to 2.0 as the disorder is weakened. The size of the largest burst increases as 
a power law of the strength of disorder with an exponent 2/3 and saturates for strongly disordered 
bundles. 



PACS numbers: 62.20.Mk, 81.40.Np, 05.90.+m, 02.50.-r 



I. INTRODUCTION 



The damage and fracture of disordered materials is 
an important scientific and technological problem which 
has attracted an intensive research during the past years 
U S 9 ■ Theoretical and experimental studies have re- 
vealed that at the beginning of the loading process of 
highly disordered materials, first micro-cracks nucleate 
randomly, covering the entire volume of the specimen 
without any spatial correlations [l|, 0, H, 0] . Approaching 
the critical load, localization occurs resulting in a single 
growing crack along which the specimen falls apart. In a 
composite system of two solid blocks glued together along 
an interface, the damage usually concentrates along the 
weak plane of the glue |f| H, 0, 0] . Loading such compos- 
ites, interface crack propagation occurs which has also 
been found to be a complex sequence of crack growth 
and arrest with interestin g sp atial and temporal fluctu- 
ations [1, 0, 0, 0, EH fill Il2j. The crackling noise ac- 
companying the failure of disordered systems (bulk or 
interface cracking) can be recorded in the form of a com- 
plicated trail of signals whose analysis provides important 
information about the microscopic dynamics of damaging 
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Fiber bundle models (FBM) are one of the most im- 
portant theoretical approaches to the progressive damage 
of disordered materials. During the last decade FBMs 
have provided a deep insight into the collective nature 
of the microscopic dynamics and statistical properties of 
degradation phenomena. Recently, FBMs have also been 
applied to study interfacial failure of glued solid blocks 
under shear loading 0, 0] and wear |(| 0, [Hj]. Inter- 
esting novel results have been obtained on the temporal 
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and spatial fluctuations of local breakings which precede 
macroscopic failure, and on the analogy of fracture with 
phase transitions and critical phenomena 0, 0, 0, 0, Ell • 
In this paper we study the damage process of a fiber 
bundle in a wedge-shape geometry, which provides a con- 
stant strain gradient of fibers. To obtain a simple rep- 
resentation of the geometry and loading conditions, we 
consider a bar subject to three-point bending. The bar 
is modelled as two rigid blocks coupled together by an 
elastic interface which is then discretized by a bundle 
of parallel fibers. Deformation and damage of the bar 
is concentrated in the interface resulting in a linear de- 
formation profile of fibers, while the two blocks remain 
intact. Besides interfacial failure, the model provides the 
mean field limit of the failure of disordered materials un- 
der three-point bending. Varying the amount of disorder 
of fibers, we can control the strength of non-linearity be- 
fore macroscopic failure, and hence, the type of fracture 
(brittle-quasi-brittle) of the bundle flElmlH}. We focus 
on the progressive damage of the fiber bundle analyzing 
the damage profile, crack formation, and bursts of local 
breakings. We find that for weak disorder an interface 
crack is formed by a continuous region of failed fibers. 
Ahead the interface crack, a process zone develops which 
proved to shrink with increasing deformation making the 
crack tip sharper. For strong disorder the failure of the 
bundle occurs due to a spatially random sequence of local 
breakings. Very interestingly we find that the size distri- 
bution of bursts is a power law whose exponent shows a 
crossover from exponent 2.5 to 2.0 when the strength of 
disorder is lowered. 

We demonstrate that the novel results of our model 
calculations are the consequence of the strain gradient; in 
the homogeneous case of zero gradient our model recovers 
all recent results of FBMs with varying threshold disorder 
[H, [H, [H [M [U HI] . Our results imply that interfacial 
fracture problems can lead to novel universality classes 
of breakdown phenomena. 
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of fibers U can be expressed as a function of 8 as 




FIG. 1: The geometrical layout of the model. The two rigid 
blocks of side length a and b are glued together by an interface 
of width lo which is discretized in terms of elastic fibers. The 
specimen suffers deflection 8 under the action of the external 
force F exerted in the middle of the bar. The wedge shaped 
opening of the interface results in a linear deformation profile 
of fibers. 



II. THE MODEL 



In order to obtain a fiber bundle with a linear defor- 
mation profile, we construct a simple model for the load- 
ing of an elastic bar of rectangular shape by an exter- 
nal force exerted perpendicular to the longer side of the 
bar in the middle. For simplicity, in the model the bar 
is composed of two rigid blocks of side lengths a and b 
which are glued together by an elastic interface of width 
lo, where Iq <C b holds, see Fig. [T] The interface re- 
gion can deform and suffer breaking under deflection of 
the specimen while the two rigid blocks remain intact. 
Bending of the specimen is performed such that the two 
blocks undergo rigid rotation about their outer upper cor- 
ner concentrating the deformation in the interface layer. 
We discretise the interface in terms of elastic fibers of 
number N and length Iq which are placed equidistantly 
between the two blocks. The fibers do not have bending 
rigidity, they can undergo only stretching deformation 
characterized by the same value of the Young modulus 
E. During the bending of the specimen, the fibers can 
support only a finite deformation, i.e. if the local de- 
formation £j of fiber i exceeds a threshold value e\ the 
fiber breaks and a micro-crack nucleates in the interface. 
The disordered properties of the material are represented 
by the randomness of the breaking thresholds which 
are independent identically distributed random variables 
with a probability density p{e) and cumulative distribu- 
tion P{e c ) = f„ E p(x)dx. The rigidity of the two rotat- 
ing blocks implies that the macroscopic deformation of 
the specimen can be characterized by a single variable 
8 which denotes the deflection of the middle of the bar 
from the original position, see Fig. [TJ 

It can be seen in Fig.[T]that under bending of the spec- 
imen the interface opens resulting in an increasing elon- 
gation of fibers from top to bottom. The actual length 
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where lx = l + 2 (b — \Jb 2 — J 2 ) is the length of fiber 
index i = 1 at the top of the bar. It follows that also the 
elongation AU and longitudinal strain Si of fibers increase 
linearly as a function of their position i 
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Equilibrium of the system is obtained when the total mo- 
mentum of forces with respect to the clamping points is 
zero. During the deformation process, those fibers which 
exceed their threshold value break, i.e. they are removed 
from the interface. Since I — P(£,(<5)) is the probability 
that the interface element of index i remained intact un- 
der the externally imposed deformation 8, based on the 
equilibrium condition, the constitutive equation a(8) of 
the deflected bar can be cast in the form 
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x [l-P( £i (S))]Eei(S), 



where L = 2b + lo is the overall length of the bar and the 
sum goes over all the fibers. On the right hand side £i(5) 
should be substituted from Eqs. (|2I3[) . The above equa- 
tions describe the macroscopic response of a fiber bundle 
which has a linear deformation profile. In the following, 
for the explicit calculations the geometrical parameters 
were set as a = 1 , b = 2.5, and Iq = 0.1. 

The amount of disorder of the failure thresholds e\ has 
a substantial effect on the macroscopic response of the 
fiber bundle In the limiting case of zero disorder, 

i.e. when all the fibers have the same breaking threshold 
el = Eq, the failure of the bundle starts at the bottom 
of the interface where the stretching deformation is the 
highest and proceeds upwards as 8 is increased. It can be 
seen in Fig [2] that the corresponding constitutive curve 
is sharply peaked. The critical deformation 8 C defined by 
the peak position corresponds to the instant of the first 
fiber breaking Eq = £jv(5 c ). Beyond the peak stress, a 
rapidly decreases due to the gradual breaking of fibers as 
8 increases. In order to study how the behavior of the 
system changes when the amount of disorder of fibers is 
varied, we consider a uniform distribution for the break- 
ing thresholds over an interval of width 2A centered at 
the value eg 



P(e c ) = 



e c - (e c - A) 
2A 



A < e c < el 



A. (5) 



The strength of disorder of the breaking thresholds is 
characterized by A, while eg sets the scale of fiber 
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FIG. 2: The constitutive curve a(S) of the interface com- 
posed of N = 10000 fibers varying the width of the threshold 
distribution A at the fixed value of Eg = 0.01. For A = the 
constitutive curve a(S) has a sharp peak which gets rounded 
and develops into a quadratic maximum when A is increased. 



strength. The width A can be varied over the interval 
[0, eg] where the limits A = and A = eg corresponds 
to zero disorder and the strongest disorder, respectively. 
Fig. [2] shows that increasing A, the peak of the consti- 
tutive curve gets more and more rounded and develops 
into a quadratic maximum. The maximum of a(S) is 
preceded by a longer and longer non-linear regime due 
to the breaking of fibers, so that for A — ► eg the linear 
behavior prevails only for small deformations 6 — > 0. On 
the microlevel this process is accompanied by the ran- 
domization of the breaking sequence of fibers along the 
interface, i.e. for A^0 fibers do not simply break in the 
decreasing order of their index i (from bottom to top of 
the bar). 



III. SIMULATION TECHNIQUES 

The complete constitutive curve of the system pre- 
sented in Fig. [2] can only be recovered by deformation 
controlled loading. When 5 is controlled externally, the 
local load on the fibers is solely determined by the exter- 
nally imposed deformation, so that there is no load re- 
distribution after fiber breaking. Under stress controlled 
conditions, the breaking of fibers is followed by the re- 
distribution of load over the intact ones. Due to the 
wedge shape of the deformed interface, at a given exter- 
nal load a, the load on the fibers linearly increases from 
top to bottom. It has the consequence that in spite of 
the rigidity of the two solid blocks, the load redistribu- 
tion following fiber failure differs from the usual equal 
load sharing (ELS) approximation commonly used for 



the study of parallel bundles of fibers @, H [U El E3 • 
The rigid surfaces, however, ensure that the load is re- 
distributed globally in such a way that the excess load 
received by an intact fiber depends on its position along 
the interface but not on its distance from the failed one. 
This implies that no stress enhancement arises in the 
vicinity of the failed fibers as in the case of local load 
sharing approximation (LLS) of fiber bundles (2l.l2ll.l24|. 
Our fiber bundle model provides the mean field limit of 
the damage and fracture of disordered materials under 
three-point bending conditions and also represents an in- 
teresting interface rupture problem. 

In order to analyze the microscopic damage mechanism 
of FBMs with a constant strain gradient, we worked out 
an efficient simulation technique for a sample where the 
interface is composed of N fibers with breaking thresh- 
olds i = 1, . . . ,N sampled from the probability distri- 
bution Eq. ([5]). Substituting the breaking thresholds e\ 
on the left hand side of Eq. ([2]) and inverting it for S, we 
can determine the value of the macroscopic deformation 
parameter 5^ — S(e?, i), i = 1, . . . , N at which the fibers 
break. Of course, d% is a function of both the position 
of the fiber i along the interface and the local breaking 
threshold s\. During the loading process the fibers break 
in the increasing order of their critical macroscopic de- 
formation 6f, which can be a randomized sequence of the 
fibers' position i. The computer simulation of the loading 
process proceeds as follows: after generating the breaking 
thresholds of fibers el we determine the corresponding 
critical deflections 5° and sort them into increasing or- 
der. The constitutive curve of the sample can be simply 
obtained by calculating the load needed to achieve the 
deformation Sf after the breaking of the first i — 1 fiber 
remaining only Ni nta ct — N — (i — 1) intact elements. 
Between the breaking of the (i — l)th and the ith fibers, 
the constitutive equation of the system takes the form 
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where the prime indicates that the summation is re- 
stricted to indices of intact fibers (which are not nec- 
essarily consecutive integers). Note that in Eq. ([6]) the 
value of S falls in the range 8^_ x < 5 < 6f. 

Performing stress controlled experiments, after the 
breaking of a fiber the deformation of the specimen can 
freely change resulting in a redistribution of load over 
the intact fibers. The excess load taken up by the intact 
fibers can give rise to further fiber failures which may 
trigger an entire avalanche of breakings. This avalanche 
either stops and the bar becomes stable under the exter- 
nally imposed load, or it spans the entire interface and 
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FIG. 3: Damage profile, i.e. the breaking probability Pt of 
fibers along the interface at five different values of the deflec- 
tion 5i < S 2 < 8 :i < 84, < 85 with A/eg = 0.2 and b/a = 2.5. 
The deflection Si falls in regime (^4), 82, 83,84, are in (B), and 
85 belongs to (C). The increase of the slope of the straight 
lines indicates the sharpening of the profile. The shadowed 
area highlights the process zone for 83. 



the specimen breaks into two pieces (the entire bundle 
ruptures) . In order to study numerically this microscopic 
breaking process, in the simulations first we increase the 
deformation i5 such that a single fiber breaks, i.e. S = 5 C { 
with index i\. Then the load needed to maintain this 
deformation a is calculated from Eq. ^ for A^„ tact = N 
fibers. After the breaking of fiber i\ its load has to be 
redistributed over the remaining N — 1 fibers. In order 
to determine the load of intact fibers after the removal of 
the broken one, we remove fiber i\ on the right hand side 
of Eq. © and invert the equation for 6(a) keeping the 
load a fixed. The fibers with threshold values 51 < 5(<r) 
break as a consequence of load redistribution. This it- 
eration has to be repeated under a fixed external load a 
until the breaking sequence stops or all the fibers break 
resulting in a macroscopic failure of the system. 



IV. SPATIAL EVOLUTION OF DAMAGE 

In a bundle of fibers loaded between two parallel rigid 
plates, due to the equal load sharing after fiber breaking, 
the failure of a fiber is solely determined by its break- 
ing threshold [H, 0, [H, [25j. Hence, fibers break in a 
completely random sequence without any spatial correla- 
tions. In our system, however, during the loading process 
the fibers break in the increasing order of their critical 
macroscopic deformation Sf {el , i) which depends both on 
the local breaking thresholds e\ and on the spatial posi- 
tion i of fibers. In the limiting case of zero disorder, i.e. 
A = and el — e^ the critical deformation 5f(eQ,i) is 
a monotonically decreasing function of i, which implies 
that the fibers break one-by-one starting from the bot- 
tom. This breaking sequence can be conceived as a crack 
is generated which penetrates the interface upward such 
that below the crack tip all the fibers are broken while 
above it the fibers are intact. Under strain controlled 
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FIG. 4: (Color online) Width of the process zone W/N as 
a function of the deflection 8 of the bar for several values of 
A. Computer simulations are in a perfect agreement with the 
analytic predictions Eqs. (|7I8I9[1 represented by the dots for 
the specific case of A/eq = 0.2. 



loading stable crack propagation is obtained gradually 
breaking the fibers by the strain increments. Controlling 
the external load, however, the onset of crack propaga- 
tion occurs in an unstable manner resulting in immediate 
catastrophic failure when the maximum of cr(d) is reached 
(see Fig.©. 

Increasing the strength of disorder A, the breaking se- 
quence of fibers determined by Sf becomes spatially ran- 
domized. At a given deformation <5, the fibers with Sf < 5 
have already failed. If the disorder is not too strong, 
an interesting spatial distribution of these broken fibers 
emerges: starting from the bottom of the interface a con- 
tinuous region of failed fibers develops forming a crack. 
On the opposite side, starting from the top of the inter- 
face a continuous region of intact, elongated fibers can 
be observed. The two regimes are separated by a process 
zone, which is a sparse sequence of intact and broken 
elements. To illustrate this feature, in Fig. [3] we show 
the probability Pb(i) that the fibers are broken along the 
interface for several different values of S using the thresh- 
old distribution Eq. ([5]). The probability Pb{i) that fiber 
i is broken at the deflection 8 can be obtained directly 
from the threshold distribution Pb(i) = P(ei(5)). The 
process zone is defined as the regime where for the prob- 
ability P\, of fiber breaking < Pb < 1 holds. It can be 
observed in Fig. that the process zone sharpens, i.e. 
its width decreases as the deformation S increases which 
makes the crack tip sharper as the crack advances. It 
can be obtained analytically that at a given deflection 5 
the width W of the process zone depends both on the 
strength of disorder A and on the geometrical extensions 
a, b of the specimen. For the explicit calculations, it is 
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worth considering separately the following three regimes 
of the damage profile: (A) for the breaking probability at 
the bottom of the interface it holds that P&(i = N) < 1, 
i.e. the crack has not yet developed (Si in Fig. [3]); (B) 
The process zone is completely contained by the inter- 
face P b (i = 1) = and P b (i = N) = 1 (S 2 ,S 3 ,S 4 in Fig. 
|3J), and (C) there is no intact region Pb(i = 1) > (S5 in 
Fig. [3]). The width W of the process zone can be obtained 
analytically as a function of 5 for the three cases 



(A) W = N 



e c Q -A-2b + 2Vb 2 -S 2 b 
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(B) W = --(N-l), 
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(C) W = N 



eg + A - 2b + 2^/b 2 -S 2 b 
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, (7) 

(8) 
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We also determined the width of the process zone W nu- 
merically for a system of 10 6 fibers, which is presented in 
Fig. 2] together with the corresponding analytic results. 
In Fig. 0]the curves of W(S) are composed of three dis- 
tinct parts corresponding to the regimes (A), (B), and 
(C) of Eqs. (|7l8l9p . It can be seen that for smaller values 
of A/eg first W increases and reaches a maximum where 
the crack occurs. As the crack advances, the width of the 
process zone decreases according to Eq. ((5J and finally, 
as the tip of the process zone reaches the top of the in- 
terface, W rapidly decreases as given by Eq. It is 
interesting to note that for large A values no crack can 
be identified, i.e. for the parameter set used in Fig. 3] the 
damage profile spans the entire interface max[W/./V] = 1 
when A/eg exceeds 0.46. It follows from the above ar- 
guments that the disorder of the interface can be consid- 
ered strong if the damage proceeds as a spatially random 
sequence of local breakings without the formation of a 
propagating crack. No crack can develop if at the defor- 
mation S where the top of the interface may already be 
damaged £i(S) > eg — A the bottom of the interface may 
still be intact En(S) < eg + A. Making use of Eq. ||2J| and 
assuming b 3> eg — A, the condition of strong disorder 
can be formulated as 



b A2 

a 1 



A>eg, 



(10) 



which implies that the average fiber strength eg, the 
width of the distribution A and the geometrical layout 
a, b of the specimen together determine the relevance of 
disorder. At a given value of a, b, and eg, the crossover 
point A* between weak and strong disorder can be ob- 
tained as 
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so that for A > A* no crack is formed, while for A < A* 
crack propagation occurs with a shrinking process zone 
ahead the crack tip. For the parameter values of Fig. 0] 
the crossover point A* = 0.4633 was obtained from Eq. 
(jlip in an excellent agreement with the numerical results. 
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FIG. 5: (Color online) a) Distribution of burst sizes D(s) 
varying the value of A in a broad range, b) Distributions 
for the limiting case of very weak disorder. Simulations were 
carried out with N = 10 6 fibers averaging over 10 3 samples. 
The geometrical layout of the sample was b/a — 2.5, and the 
scale parameter of fibers' strength had the value Eq = 0.01. 
Straight lines of slope 2.5 and 2.0 are drawn as reference. 



V. BURSTS OF FIBER BREAKINGS 

In order to characterize the damage process of the fiber 
bundle under stress controlled conditions, we determined 
the distribution of burst sizes s of fiber breakings varying 
the width A of the disorder distribution. Simulations 
were carried out by increasing the external load to break 
a single fiber and following the cascading fiber breakings 
with the algorithm discussed in Sec. IIIII The avalanche 
size distributions D(s) are presented in Fig. O varying the 
amount of disorder in a broad range. It can be observed 
in Fig. [5] that D(s) has a power law form 



D(s) 



(12) 



at any finite value of A with an exponential cutoff at large 
avalanches. Simulations revealed an interesting change of 
the value of the exponent r as the strength of disorder is 
varied: for strong disorder the value of r coincides with 
the mean field (equal load sharing) ex pon ent of the clas- 
sical parallel bundle of fibers r ps 2.5 jl3l . [26[. However, 
as the disorder is weakened, the distribution exhibits a 
crossover to another power law with a lower exponent 
r ps 2.0. To better illustrate this effect, in Fig.[5j>) burst 
distributions are shown separately for the limiting case 
of very weak disorder. The numerical results are well 
described by a power law with an exponent 2.0. 

It is interesting to note that the size of the largest burst 
Smax has a strong dependence on the value of A (see 
Fig. EU&)), namely, as the strength of disorder is reduced 
the largest avalanche decreases. To obtain a quantitative 
characterization of this effect, Fig. [6] presents the average 
size of the largest bursts (s max ) as a function of A, where 
a power law dependence is evidenced 



A Q 



(13) 



for the case of weak disorder A < A*. As A exceeds A*, 
the largest avalanche s max reaches a maximum and levels 
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FIG. 7: (Color online) a) Burst size distributions for the lim- 
iting case of zero strain gradient (a — 0.0). There exists a 
finite critical disorder A c ~ 0.258, where a crossover occurs 
from the exponent 5/2 to a lower value 3/2. (b) At any finite 
value of a > criticality is obtained only for A — > and the 
crossover exponent becomes larger. 



FIG. 6: Average size of the largest avalanche as a function 
of A for several values of eg. A straight line of slope 2/3 is 
drawn to guide the eye. 



off (see Fig. [5]). Computer simulations revealed that the 
exponent a has a universal value a = 2/3 indicated by 
the straight line drawn in Fig. [5] to guide the eye. The 
value of A where (s max } saturates is in a good agreement 
with the corresponding value of A* estimated from Eq. 

cm. 

It is important to emphasize that in our system no crit- 
ical disorder distribution can be identified in the sense 
defined in Refs. 0, M, M, HI H, HI . For equal load 
sharing fiber bundles the threshold distribution is consid- 
ered to be critical if the breaking of the weakest fiber gives 
rise to an immediate macroscopic failure of the system. 
For uniformly distributed threshold values in the inter- 
Xg, 1] the distribution becomes critical for xq — > 0.5 
lli^SHIl^l. It has been pointed out that ap- 
proaching the critical disorder, the macroscopic response 
of the system becomes perfectly brittle and on the micro- 
level the size distribution of bursts exhibits a crossover 
from a power law of exponent 2.5 to a significantly lower 
one 1.5 [HI, [2l|. Using our terminology, such critical 
behavior in a bundle of fibers loaded between two par- 
allel rigid plates, i.e. in the case of equal load sharing, 
should be obtained for A — > Sq/3, however, computer 
simulations revealed a sudden collapse of our FBM with 
constant strain gradient solely at A = 0. 

The reason for the missing critical state in our model 
is the inhomogeneity of the load of fibers along the inter- 
face. Eq. ([2]) shows that at any deflection 5 the load of 
intact fibers linearly increases from top to bottom of the 
interface where the gradient, i.e. the strength of inhomo- 
geneity, is determined by the geometry of the system a/b. 
Setting the cross section of the specimen a to zero a = 0, 
the positional dependence of £j disappears in Eqs. 



and formally all the fibers keep the same load determined 



by 5. We carried out computer simulations in the limit- 
ing case a = varying the strength of disorder A. We 
indeed find that in this homogeneous case a critical state 
arises at A c /eg ss 0.258 such that for A < A c a single 
fiber breaking triggers a catastrophic avalanche. The cor- 
responding numerical results are presented in Fig. [7](a), 
where it can be seen that approaching A c from above, 
the burst size distribution exhibits a crossover from the 
exponent r = 5/2 to a lower value r = 3/2 in agreement 
with the predictions of Refs. [HI, H3, [HI . The deviation 
of A c from £q/3 arises due to the non- linear terms of 5 
in Eqs. (|2I4|) . It can also be observed in the figure that 
in the homogeneous case the largest burst s max increases 
as criticality is approached in agreement with Ref. [2l|. 

However, at any finite value of a, i.e. in the presence of 
a strain gradient, the picture drastically changes: critical- 
ity occurs solely in the limit A — > so that catastrophic 
failure is always preceded by avalanches with a power law 
distribution but the cutoff avalanche size goes to zero as 
a power law of A when A is decreased. Fig. [7f)) presents 
simulation results obtained with 10 6 fibers at the value 
a = 0.01, which results in a much lower strain gradient 
than in Fig. [5] Comparing the burst size distributions of 
different values of A for the homogeneous (Fig. Wp)) and 
inhomogeneous (Fig.©)) cases, a clear difference can be 
observed. The higher value of the crossover exponent 
t = 2 compared to t = 3/2 of the homogeneous case 
shows that in the absence of a finite critical disorder, the 
large avalanches are less dominating in the distributions 
D(s). 



VI. DISCUSSION 

We presented a detailed study of the progressive dam- 
age and fracture of fiber bundles in a wedge-shape geom- 
etry which provides a linear deformation profile for fibers. 
For a simple representation of the geometrical and load- 
ing conditions of the system, we considered a bar sub- 
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jected to three point bending. The bar is composed of 
two rigid blocks coupled by an elastic interface which is 
then discretized in terms of parallel fibers. We showed 
that in the limit of zero disorder of fibers' strength the 
bundle has a perfectly brittle macroscopic response, i.e. 
under stress controlled loading global failure occurs as a 
sudden collapse of the system without any precursory ac- 
tivity, furthermore, fibers break in a completely ordered 
sequence from bottom to top of the wedge creating an 
instable crack with a sharp tip. The relevance of disor- 
der is determined together by the parameters (mean and 
width) of the strength distribution of fibers and by the 
geometrical layout of the wedge. We demonstrated that a 
propagating interface crack can only be defined for weak 
disorder. Ahead the crack a process zone is formed whose 
width decreases with increasing deformation making the 
crack tip sharper as the crack advances. For strong disor- 
der a spatially random sequence of local breakings occur 
along the entire bundle. 

The breaking of single fibers can trigger cascades of 
breaking events. The size distribution of these bursts is 
found to be a power law with an interesting crossover 
effect as the strength of disorder is varied: for strong 
disorder the mean field exponent r = 5/2 of equal load 
sharing fiber bundles is recovered indicating the complete 
randomness of the failure process. However, for weak dis- 
order where a propagating crack with a process zone de- 
velops, a lower exponent r = 2.0 is obtained. In the weak 
disorder regime the largest burst increases as a power law 
of the width of the disorder distribution with an expo- 
nent a = 2/3. We showed that in the limit of zero strain 
gradient our calculations reproduce the crossover of burst 
exponents from 5/2 to 3/2 predicted recently. The novel 
features of our system originate from the finite strain 



gradient, i.e. from the geometrical constraints of fibers 
which naturally occur, for instance, at interfacial frac- 
ture problems. Our model provides also the mean field 
limit of the damage and fracture of disordered materials 
under three-point bending. The results imply that the 
statistical physics of interfacial rupture can reveal novel 
universality classes of breakdown phenomena. 

Recently, it was found experimentally that the crack- 
ing of a bar under three point bending proceeds in bursts 
which are characterized by power law distributions [3lj ] - 
The experiments showed that the exponents of the ampli- 
tude, area and energy distribution of magnetic emission 
signals recorded during the fracture process of ferromag- 
netic materials are sensitive to the type of fracture, i.e. 
the noise spectra of ductile materials are characterized 
by higher exponents than the brittle ones. The bound- 
ary and loading conditions ensured in the experiments 
that the damage localizes to a relatively thin layer of the 
specimen giving rise to a single growing crack so that the 
crackling noise measured during the loading process char- 
acterizes the crack propagation. Note that in our model 
the crossover to a lower exponent of burst sizes with de- 
creasing disorder is accompanied on the macro-level by 
an increasing degree of brittleness showing that this sim- 
ple mean field approach can qualitatively account for the 
changing properties of crackling noise observed experi- 
mentally [31] . 
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